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The notion of zitterbewegung is a long-standing prediction of relativistic quantum mechanics. Here 
we extend earlier theoretical studies on this phenomenon for the case of III-V zinc-blende semicon- 
ductors which exhibit particularly strong spin-orbit coupling. This property makes nanostructures 
made of these materials very favorable systems for possible experimental observations of zitterbe- 
wegung. Our investigations include electrons in n-doped quantum wells under the influence of both 
Rashba and Dresselhaus spin-orbit interaction, and also the two-dimensional hole gas. Moreover, we 
give a detailed anaysis of electron zitterbewegung in quantum wires which appear to be particularly 
suited for experimentally observing this effect. 

PACS numbers: 



I. INTRODUCTION 

The spin degree of freedom of electrons in semicon- 
ductor nanostructures is one of the central subjects in 
the growing field of spin electronics 0. The latter key 
word describes the whole variety of efforts and propos- 
als for using the electron spin instead, or in combination 
with, its charge for information processing, or, even more 
ambitious, quantum information processing. A central is- 
sue is the possibility of electrical control electron spins, 
which avoids many difficulties arising from applying and 
gating magnetic fields. Zinc-blende III-V semiconductors 
show a particularly strong spin-orbit interaction and are 
therefore natural candidates with respect to the above 
goals. 

In the present paper we extend earlier results on the 
relativistic effect of zitterbewegung in systems of the 
above type @|. In this context we also discuss related 
recent work by other authors 0, 0, IE IE 111 ■ Our in- 
vestigations include electrons in n-doped quantum wells 
under the influence of both Rashba and Dresselhaus spin- 
orbit interaction, and holes in two-dimensional geometry. 
Phenomena related to spin-orbit coupling in such systems 
are presently studied very intensively also in the context 
of the intrinsic spin- Hall effect ^LflUJ- Particular at- 
tention is paid to the case of electron zitterbewegung in 
quantum wires which appear to be good candidates for 
experimental investigations of this effect. 

The coupling between the orbital and the spin degree 
of freedom of electrons is a relativistic effect described 
by the Dirac equation and its nonrelativistic expansion 
in powers of the inverse speed of light c 12] . In second 
order one obtains, apart from two spin-independent con- 
tributions, the following well-known spin-orbit coupling 
term, 

H so = -^s- (w x -P-V (1) 

2moc z \ mo ) 

where mo is the bare mass of the electron, s, p its spin and 



momentum, respectively, and V is some applied external 
potential. On the other hand, the free Dirac equation, 
V = 0, has two dispersion branches with positive and 
negative energy, 

e(p) = ±y/m? c* + c?p>, (2) 

which are separated by an energy gap of 2moc 2 « lMeV. 
In particular, the nonrelativistic expansion of the Dirac 
equation quoted above can be seen as a method of sys- 
tematically including the effects of the negative-energy 
solutions on the states of pos itive energy starting from 
their nonrelativistic limit [l^- Moreover, the large en- 
ergy gap 2moc 2 appears in the denominator of the right 
hand side of Eq. JQ), suppressing the effects of spin-orbit 
coupling for weakly bound electrons. 

On the other hand, the band structure of zinc-blende 
III-V semiconductors shows many formal similarities to 
the situation of free relativistic electrons, while the rele- 
vant energy scales are grossly different [IE [Tj, [IE] • For 
not too large doping of such semiconductors, one can 
concentrate on the band structure around the T point. 
Here one has a parabolic s-type conduction band and a 
p-type valence band consisting of the well-known disper- 
sion branches for heavy and light holes, and the split-off 
band. However, the gap between conduction and valence 
band is of order leV or smaller. This heuristic argument 
makes plausible that spin-orbit coupling is an important 
effect in III-V semiconductors which actually lies at the 
very heart of the field of semiconductor spintronics. 

Dating back to a seminal paper by Schrodinger |T3.IT7| 
from 1930, the notion of zitterbewegung has been a long- 
standing theoretical prediction of relativistic quantum 
mechanics. In the free Dirac equation, this oscillatory 
quantum motion occurs for particle wave packets which 
are superpositions of both solutions of positive and neg- 
ative energy. Thus, the dominant frequency of this dy- 
namics is given by 2moC 2 /h which is of order 10 20 Hz. 
Moreover, the length scale of this motion is given by the 
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Compton wave length h/m^c of the free electron. There- 
fore, in order to experimentally observe the zitterbewe- 
gung of free electrons one would need to confine these 
objects on a length scale of a few picometers. Now it fol- 
lows form general uncertainty arguments that such a spa- 
tial confinement leads to an energy scale where electron- 
positron pair production plays a serious and detrimental 
role. These arguments have led many authors to the 
opinion that the zitterbewegung of electrons is impossible 
to observe, see e.g. Ref. 18]. 

Most recently, we have theoretically investigated zit- 
terbewegung in III-V zinc-blende semiconductors and de- 
veloped a proposal for its experimental observation. As 
already mentioned, the mathematical treatment of effec- 
tive band structure models relevant to such systems on 
the one hand, and the nonrelativistic expansion of the 
Dirac equation on the other hand show many formal sim- 
ilarities 0,0,01- In fact, under certain aspects, the 
low-energy band structure of such semiconductors around 
the Fermi level can be viewed as a model for truly rela- 
tivistic electrons, but with energy and length scales which 
are much more favourable for observing effects like zit- 
terbewegung. 

This paper is organized as follows. In section [H] we 
study zitterbewegung in the two : dimensional electron gas. 
We extend previous results Q to the situation where 
spin-orbit coupling of both the Rashba and Dresselhaus 
type is present. Moreover, we give a deeper discussion 
of zitterbewegung in quantum wires and our related ex- 
perimental proposal. Section IlIII is devoted to the same 
phenomenon for heavy holes in two-dimensional quantum 
wells. We close with conclusions in section llVl 
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II. ZITTERBEWEGUNG IN THE 
TWO-DIMENSIONAL ELECTRON GAS AND IN 
QUANTUM WIRES 

A. The two-dimensional electron gas 

For conduction band electrons in n-doped quantum 
wells of zinc-blende semiconductor structures the dom- 
inant effects of spin-orbit interaction can be described in 
terms of two effective contributions to the Hamiltonian. 
On of them is the Rashba spin-orbit term which is 
due to the inversion-asymmetry of the confining potential 
and has the form 



n R = - ( Px <jy 



Py° X ) , 



(3) 



where p is the momentum of the electron confined in 
a two-dimensional geometry, and a the vector of Pauli 
matrices. The coefficient a is tunable in strength by 
the external gate perpendicular to the plane of the two- 
dimensional electron gas. The other contribution is the 
Dresselhaus spin-orbit term which is present in semicon- 
ductors lacking bulk inversion symmetry 20]. When re- 
stricted to a two-dimensional semiconductor nanostruc- 
ture grown along the [001] direction this coupling is of 



(4) 



where the coefficient (5 is determined by the semiconduc- 
tor material and the geometry of the sample. These two 
contributions to the effective Hamiltonian have also an 
interesting interplay |U 0] • 

We now consider the single-particle Hamiltonian of a 
free electron under the influence of spin-orbit coupling of 
both the Rashba and the Dresselhaus type, 



->-> 

H = 1- Hr + Hd , 

Zm 



(5) 



where m is the effective band mass. The components of 
the time-dependent position operator 

r H (t) = e mt/h f{0)e- tHt/h (6) 

in the Heisenberg picture read explicitly 

x H (t) = x{Q) + + ( %y - ^a x )t 
m \n a J 

a 2 -/3 2 ( ( 2 A \\ h z 

+ ^^(l A *- sin (l A< )) s ' (7) 

ya (t) = x {o) + p -it-(^-Ly)t 

- ^.(|a<-™(|a«)) S (8) 



with 



and 



A 2 (p; a, 13) = (a 2 + (3 2 ) p 2 + 4a(3p xPy 



(9) 



E(p,a;a,/3) = (p x a x + p y a v ) 

+ P( Px <jy +Py<T X )) . (10) 



Here the operators p and a on the right hand sides are in 
the Schrodinger picture and therefore time-independent. 

The oscillatory terms on the right hand sides of 
Eqs. J7J, JHJl can be viewed as the zitterbewegung the elec- 
tron performs under the influence of spin-orbit coupling. 
This oscillatory quantum motion vanishes if relativistic 
effects are absent, a — f3 — 0. The contributions linear 
in time t in the first lines in Eqs. 0,® ar e just propor- 
tional to the velocity 



v = -[H,r] 



(11) 
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which is in the presence of spin-orbit coupling a spin- 
dependent operator. The expressions Qj© contain the 
results given in Ref. 2] for pure Rashba or Dresselhaus 
coupling as special cases. For, instance, if only Rashba 
coupling is present (f3 = 0), one can express the position 
operator in the Heisenberg picture as 



x H (t) = x (0) + ^t+^^a z [ 1-cos 

777 p z 2 1 



x \ ( 2a P. ■ ( 2a P.\\ 



y H (t) 



(12) 



(13) 



The case a = ±/3 is particular |23l |24| . As seen from 
Eqs. Ql,JHJ, the oscillatory part of the time-dependent 
position operator, i.e. the zitterbewegung vanishes due 
to the prefactor (a 2 — (3 2 ). Physically, this observation 
results from an additional conserved quantity arising at 
this point pll23|. 

It is straightforward to evaluate the above time- 
dependent position operators within Gaussian wave pack- 
ets. For simplicity we concentrate on the case of pure 
Rashba coupling described by Eqs We con- 

sider a Gaussian wave packet with initial spin polariza- 
tion along the z-direction perpendicular to the quantum 
well, 



1 d 
2^~~ 



TT 



(14) 



Clearly we have (ip\r\tp) = 0, {ip\0[ip) = hko, and the 
variances of the position and momentum operators are 
(Ax) 2 = (Ay) 2 = d 2 /2, (A Px f = (A Py f = h 2 /2d 2 . 
Thus, the group velocity of the wave packet is given by 
hko/m , while its spatial width is described by the pa- 
rameter d with the minimum uncertainty product typical 
for Gaussian wave packets, Ap x Ax — Ap y Ay = h/2. 
A direct calculation gives 



hk 0x 



d 

2^ 



-rflk 2 



271 



dip sin if 



171 -n Jq 

q 2 +2dq(kox cos (p+ko v sin ip) 



dq e 



(i 



/2aq 

1 — cos 1 

\ nd 



(15) 



In the above expression q is a dimensionless integration 
variable. The remaining integration over the polar angle 
<p gives a vanishing result if ko y = 0, i.e. if the group 



velocity is along the x-direction. More generally, one 
finds that 



(^•Fff(t)/fco|V) = — *: 
777, 



(16) 



which means that the zitterbewegung is always perpendic- 
ular to the group velocity of the wave packet. The same 
observation is made for the case of pure Dresselhaus cou- 
pling, a = 0. Note that for pure Rashba coupling the 
sum of the z-component of the orbital angular momen- 
tum I = fx p and the z-component of the spin s — Ha/2 
is a conserved quantity, \H, l z + s z ] = 0, while for pure 
Dresselhaus coupling we have [H, l z — s z ] — 0. In the 
presence of these conservation laws, the zitterbewegung 
can also be interpreted as a consequence of spin rota- 
tion due to spin-orbit coupling. Consider for instance an 
electron moving along the y-dircction with its spin being 
initially aligned with the z-dircction. In the time evo- 
lution of the particle, the spin will then be rotated due 
spin-orbit coupling which requires, by virtue of the con- 
servation law, also a finite component l z to develop, i.e. 
the electron has to perform also a movement perpendic- 
ular to its group velocity. In the general case a ^ 7^ (3 
such an interpretation does not seem to be possible since 
a conserved quantity of the above kind does not exist. 

Let us turn back to the case of pure Rashba cou- 
pling. Without loss of generality we consider an electron 
wave packet moving along the y-direction, k$ x = 0. By 
expanding the exponential containing the trigonometric 
functions in Eq. (|15fl . one derives [25| 



1 



2h 



1 



-d'kf. 



00 



(dk 0y ) 



2(n+l) 



n\(n+ 1)! 



dqq 2 



(2aq + 



.(17) 



Thus, the amplitude of the zitterbewegung is proportional 
to the wave length Aoj, = 2t: /ko y of the electron motion 
perpendicular to it. In a semiconductor quantum well, 
this length can be of order a few ten nanometers, which is 
several orders of magnitude larger than the length scale of 
the zitterbewegung of free electrons given by the Compton 
wave length. Note also that the oscillatory zitterbewegung 
changes its sign if the translational motion is reversed. 

If the product dk 0y is not too large, dk 0y < 1, only 
low values of the summation index n in Eq. I|17|l lead 
to substantial contributions, and the Gaussian factor in 
the integrand suppresses contributions from large values 
of q. Thus, a typical scale of this integration variable 
is leading to sizable contributions is q ~ 1/V2. Thus, 
a typical time scale in the integrand is T = \/2irhd/a, 
and when averaging the zitterbewegung over times scales 
significantly larger than T, the cosine term drops giving 



{il>\x H (t)\il>) 



1 



2k, 



0y 



(l-e- d ' 2k '»y) 



(18) 
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i.e. the time-averaged guiding center of the wave packet 
is shifted perpendicular to its direction of motion. Note 
that the zitterbewegung is absent for ko y = |l8j . 

In the opposite case dko y 3> 1 the Gaussian approaches 
a 5-function. In this limit one finds (for k$ x = 0) 

m^wm = ^ i 1 cos (^H) • (19) 

Here the frequency of the zitterbewegung is f2 = 2ako y /H, 
and the guiding center of the wave packet is also 
shifted in the direction perpendicular to its group ve- 
locity. Note that Ml is the excitation energy between 
the two branches of the Rashba Hamiltonian Ti at a 
given momentum k — k$ y e y . Rashba spin-orbit cou- 
pling is particularly strong in InAs where values for 
the parameter a of a few 10 -11 eVm can be reached 
EE IS EHH El EH!!, leading to frequencies ft is the 
terahertz regime. This is much smaller than the typical 
frequency of the zitterbewegung of free electrons which 
is of order 10 20 Hz. For GaAs, the Rashba coefficient 
is typically an order of magnitude smaller than in InAs 
[33| and the Dresselhaus coupling plays a more impor- 
tant role summary, the zitterbewegung of 
electronic wave packets in semiconductor quantum wells 
as discussed above is characterized by amplitudes and 
frequencies which are by orders of magnitude larger and 
smaller, respectively, than it is the case for free electrons. 
This opens the perspective to experimentally observe the 
electron zitterbewegung in semiconductor nanostructures 
via terahertz methods, or using high-resolution scanning- 
probe microscopy ima ging techniques as developed and 
discussed in Refs. [37t l38j] The latter approach will be 
described in more detail in section III Bl 

The issue of zitterbewegung of electrons in III-V semi- 
conductors was also discussed recently by Zawadzki 
where the three-dimensional bulk case was considered |^ . 
The author starts from an 8 x 8 Kane model for conduc- 
tion and valence band being diagonal in the hole sector 
|l4i I39I ]. Specializing on particles moving along the z- 
direction in real space and neglecting the split-off band, 
the author obtains an effective Hamiltonian coupling only 
light holes and conduction band electron states. More- 
over, for an appropriate choice of basis this Hamiltonian 
matrix mimics the Hamiltonian of the free Dirac equa- 
tion, which enables to derive a zitterbewegung following 
directly Schrodinger's original approach (la. fl7| . How- 
ever, when reducing the underlying 8x8 Kane model in 
a systematic way to an effective 2x2 Hamiltonian for con- 
duction band electrons only, one obtains in the absence 
of an additional potential and magnetic fields to second 
order in the gap energy only a kinetic term involving an 
effective mass depending on band structure parameters 
[l4|. This is analogous to the situation of the free Dirac 
equation where, again in the absence of external fields, 
the lowest-order relativistic correction is an additional 
contribution to the kinetic energy which does not lead 
to zitterbewegung |12| . The zitterbewegung studied in the 
present work occurs, for the case of Rashba coupling, due 



to an external potential introducing structure-inversion 
asymmetry [HHg. For the case of Dresselhaus coupling, 
it stems from the bulk Dresselhaus coupling term which 
results from bulk-inversion asymmetry |20| and is not in- 
cluded in the 8x8 Kane model. In this sense, Zawadzki's 
result appears to be effect of higher order in the inverse 
gap energy which can be of importance in materials with 
particularly small gap such as InSb |39| . Moreover, the 
zitterbewegung as discussed in Ref. Q occurs always in 
in the direction of the group velocity of the particle wave 
packet, i.e. in the z-direction. This feature is clearly 
different from the zitterbewegung investigated here and 
might pose an obstacle against experimentally observing 
this effect. 



B. Quantum wires 

The zitterbewegung of an electron in a quantum well 
as described above is naturally accompanied by a broad- 
ening of the wave packet, where the dominant contribu- 
tion stems from the dispersive effective-mass term in the 
Hamiltonian. Such a broadening might pose an obstacle 
for experimentally detecting the zitterbewegung. How- 
ever, the broadening can be efficiently suppressed and 
limited if the electron moves along a quantum wire. In 
fact, the motion of electrons in quantum wells is generally 
under better control if additional lateral confinement is 
present. We therefore consider a harmonic quantum wire 
along the y-direction described by 

W = |- + \mu 2 x 2 + H R , (20) 
2m I 

where the frequency to parameterizes the confining poten- 
tial perpendicular to the wire 0, . It is instructive 
to write the Hamiltonian in the form 

H = H Q + H! (21) 

with 

H = hw(a+a+-\ + + ak 0y a x , (22) 

Hi = -i^__( a -o+)a*. (23) 

Here a, a + are the usual harmonic climbing operators, 
and fcoy is the component of the electron wave vector 
along the quantum wire. Due to the properties of the 
"mixing operator" Ti.± , analytical progress as before with- 
out employing further approximations does not seem to 
be possible. We therefore project the Hamiltonian 7i 
onto the lowest two orbital subbands. This approxima- 
tion is known to give very reasonable results for not too 
wide quantum wells |4lj| , and we will also compare its re- 
sults with a full numerical simulation of the above multi- 
band Hamiltonian. 

For a given k$ y the truncated Hilbert space is spanned 
by the states |0, f), 10,1), |1jT)j 1 1 , 4,) , where the arrows 
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denote the spin state with respect to the ^-direction, 
and and 1 stand for the ground state and the first 
excited state of the harmonic potential, respectively. In 
the above basis, the truncated Hamiltonian reads 



n 
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-ak 0y 
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-ak 0y 
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huje / 
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-ah 



-hu> 
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II K 0y 
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-ak y 



(24) 



(25) 
(26) 



being the subband energies in the absence of spin-orbit 
coupling, and en = mo? jf? is the energy scale of the 
Rashba coupling. When applying the transformation 



(27) 



the projected Hamiltonian and in turn the time evolution 
operator become block-diagonal, 




Ue 



\Ht TT + 



U + = 



h+(t) 
h-{t) 



(28) 



where 
h±(t) 

and 



1 cos 



• exp 




jmuj Huj 
M± = [ ±«\/ -^fi, 0, — — T ak 0y 



(30) 



Let us first consider an electron with a given momentum 
ko y along the wire and injected initially into the lowest 
subband of the confining potential with the spin point- 
ing upwards along the z-direction, i.e. the initial wave 
function \ip z ) — |0, t) for the x-direction is a Gaussian 
whose width is determined by the characteristic length 
A = y^h/muj of the harmonic confinement. Using 



UxU+ = — ( ° ^ 
V2\° x 



(31) 



one obtains for the above initial state the following time- 
dependent expectation value 



{i>z\XH{t)\lpz 



\/2 M+M- 

• sin ( j^^+t J sm (j^f 1 ^ 



(32) 



The amplitude of this oscillatory dynamics perpendic- 
ular to the wire direction becomes maximal when the 
resonance condition Hlu — ±2ako y is fulfilled. At that 
point we have fi^- = 0, and if fi± can be neglected com- 
pared to /ifj_ (which is the case for large enough ko y ) the 
amplitude of the zitterbewegung is approximately \ j \[2. 
This result from the truncated Hamiltonian is in excellent 
agreement with numerical simulations of the full multi- 
band system. In Fig. ^ we have plotted simulation re- 
sults for (ipzlxHitflipz) where the wave number along the 
wire is fixed to be ko y X = 5 and the Rashba parame- 
ter a is varied around the resonance condition. Clearly, 
the amplitude is maximum at resonance. Equivalent ob- 
servation are made if the Rashba coupling is fixed while 
the wave number ko y is varied. In Fig. |3 we have plot- 
ted the amplitude of the zitterbewegung as a function of 
Vl/u) = 2akQ y /huj for different values of the wave number 
ko y along the wire. In this range of parameters, the res- 
onance becomes narrower with increasing ko y , while its 
maximum value is rather independent of this quantity 
and remarkably well described by X/y/2. 

A qualitative explanation for this resonance can be 
given in terms of the decomposition H = Ho + T~li of 
the Hamiltonian. The zitterbewegung is induced by the 
perturbation 7ii which can act most efficiently if the un- 
perturbed energy levels of Ho are degenerate having op- 
posite spins. This is the case at |2a&ch/| = f^- 

We propose that electron zitterbewegung in semicon- 
ductor nanostructures as described above can be exper- 
imentally observed using high-resolution scanning-probe 
microscopy imaging techniques as developed and dis- 
cussed in Refs. (33, yg. As a possible setup, a tip can be 
moved along the wire and centered in its middle. For an 
appropriate biasing of the tip, the electron density at its 
location is partially depleted leading to a reduced conduc- 
tance of the wire. Since the amplitude of the zitterbewe- 
gung reflects the electron density in the center of the wire, 
the zitterbewegung will induce beatings in the wire con- 
ductance as a function of the tip position. These beatings 
are most pronounced at the resonance, see Fig. ^ Note 
that the oscillations shown there as a function of time 
can be easily converted to the real-space y-coordinate by 
multiplying the abscissa by hko y /m. Moreover, for tun- 
ing the system to the resonance condition |2afcoj,| = huj, 
at least two parameters can be varied expertimcntally: 
The group velocity of the injected electron along the wire 
given by ko y , and the Rashba parameter a which is tun- 
able by a gate voltage across the quantum well. Thus, 
quantum wires defined in InAs quantum wells are favor- 
able systems for experiments of the above kind, since this 
material can exhibit a quite large Rashba coupling but 
has a comparatively small Dresselhaus term. The group 
velocity can be varied by changing the gate to the two- 
dimensional electron gas. This alters the global electron 
density and therefore also the wave vector for motion 
along the wire. Generally we expect spin-orbit effects in 
STM experiments to be more pronounced in the presence 
of additional confinement such as in a quantum wire. 



6 



The present considerations concentrate on the case of 
pure Rashba coupling neglecting a possible Dresselhaus 
contribution. Very analogous observations as above can 
be made for pure Dresselhaus coupling, while in the case 
of both couplings being present the analytical theory be- 
comes technically more involved. We note that the Dres- 
selhaus coupling, differently from the Rashba term, can- 
not be tuned by an external gate. Therefore, materials 
with pronounced Rashba coupling are favorable for tun- 
ing the system to the above resonance condition. 

Let us now analyze the situation when the electron is 
injected into the lowest subband of the wire, but its spin 
is not aligned with the z-direction. If the spin points 
along the x-direction, no zitterbewegung occurs, 



(ip x \x H (t)\ip x ) = 0. 



(33) 



This is a property of both the full multi-band model and 
the truncated Hamiltonian and follows from symmetry 
considerations: Under a reflection in the yz-plane x, p x 
and the y- and z-component of the spin change sign while 
the other components of spin and momentum remain un- 
changed. Thus, the Hamiltonian is invariant under this 
operation, a property which is also shared by the initial 
state \ip x )- Therefore, the expectation value of xn{t) has 
to be equal to its negative and is consequently zero. Fi- 
nally, if the spin points initially along the ^-direction one 
finds from the truncated Hamiltonian 



(i/jy\x H {t)\ip v ) 



A 1 
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, 1 \ • f 1 

■ sin ( Tfi+t I sm I J.^-^ 



/*+ (1 \ . f 1 
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H+ \h J \h 

V- . fl \ ( 1 

— sm — \i+t cos —fi-t 



(34) 



Thus, the zitterbewegung also occurs if the electron spin 
is initially aligned along the wire direction. 

It is instructive to also investigate the dynamics of the 
spin degree of freedom as the electron passes along the 
wire. For a situation where the spin is pointing again 
initially in the z-direction, the truncated Hamiltonian 
leads to the following expressions: 



(^K(i)|^ 



COS" I J^fi+t 
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— ^ sm U M+ 
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,(36) 



( 1 \ 1 

cos — /i+t cos —\X-t 
\h J \H 
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• sm I Tft+t 1 sm I J.^-^ 



(37) 



Interestingly the expectation value {ip z \<j^(t)\il) z ) 
shows a particular behavior at the resonance frw — 
±2ako y . Here we have = 0, and if /x± can again 
be neglected compared to (as it is the case for the 
choice of parameters used in Figs. I1I2|| . this expectation 
value is approximately given by 



(V>»K(t)|^ 



± i 



(38) 



Thus, the time dependence is, to a very good degree of 
approximation, governed by a single frequency. This re- 
markable result is also confirmed by numerical simula- 
tions of the full multiband model shown in Fig. [3] 

Finally, if the electron spin points initially in the x- 
direction, we have 



+ 



-_H+t 



/M K39) 



while the expectation values of the other two spin com- 
ponents strictly vanish. The latter result is also true 
for the full multi-band model and follows from the same 
symmetry considerations as above. 

Let us finally summarize other recent theoretical de- 
velopments pertaining to the issue of zitterbewegung in 
quantum wires. The electron dynamics in ballistic quan- 
tum wires in the presence of spin-orbit coupling were also 
recently analyzed by Nikolic, Zarba, and Welack concen- 
trating on transverse forces on the electron induced by 
spin-orbit coupling [(| . A study similar in spirit was car- 
ried out by Shen, who also made a connection between 
transverse forces due to spin-orbit interaction and the 
occurrence of zitterbewegung jjj. Lee and Bruder stud- 
ied recently quantum wires with spin-orbit coupling of 
both the Rashba and Dresselhaus type p|. Their ap- 
proach is mainly numerical and concentrates on mapping 
out charge- and spin-density modulations in ferromagnet- 
semiconductor single-junction quantum wires. The shape 
of these modulations is explained in terms of the symme- 
try properties of the eigenstates of the wire. Finally, Za- 
wadzki has most recently analysed the band structure of 
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narrow-gap single-wall carbon nanotubes making a con- 
nections to relativistic effects in the free Dirac equation 
. Similarly to Ref. Q , the kind of zitterhewegung pre- 
dicted from these investigations is along the group ve- 
locity of the electron wave packet, i.e. along the wire, 
which is considered there as a strictly one-dimensional 
system. This is in contrast to the present study where 
we consider a quantum wire of finite width and obtain a 
zitterhewegung of electronic wave packets perpendicular 
to the wire direction. 



III. HOLES IN A TWO-DIMENSIONAL 
QUANTUM WELL 



We now turn to the case of holes in the p-type va- 
lence band of 111-V semiconductors as opposed to s-type 
conduction band electrons studied so far. We note that 
Jiang et al. have very recently performed a semiclassical 
study of the time evolution of holes in three-dimensional 
bulk systems under the influence of a homogeneous elec- 
tric field • This investigation was motivated by the re- 
cent prediction of the intrinsic spin-Hall effect 
Here we shall analyse the full quantum time evolution of 
heavy-hole states in quantum wells, a scenario for which 
spin-Hall transport was most recently predicted [Til ] and 
experimentally reported |42| . 

At low temperatures and for not too wide wells, only 
heavy-hole states are occupied with their angular mo- 
mentum pointing predominantly along the growth di- 
rection [43j, corresponding to the total angular momen- 
tum states ±3/2. Due to this constraint, the effects of 
structure-inversion asymmetry on the hole spins are tri- 
linear in the momentum, and the Hamiltonian incorpo- 
rating this type of spin-orbit coupling reads for appropri- 
ate growth directions of the quantum well 0, l44j . 



2m 2h 3 V 



p%a. 



(40) 



time-dependent position operator read 

Xa (t) = x(0) + ^i+^^Wl-cos 
m p z 2 

Px 3ft 



pry ((3PxP? - pI) <? x (Splp y - pi) ay) 



2ap 3 \ 2ap 3 
sin | . , t T^—t 



h 4 



h 4 



+ ^t(6p x JV7*+3(pg-j£)o*) , 



n 3 

y H (t) = 2/(0) 



(41) 



m p 2 2 V V h J J 



El 


m 


p 5 


~2 


( 


sin 


a 





3 \ „x 



n 3 \ n v\ 



h 3 



(42) 



Again, the zitterhewegung of a wave packet with its spin 
pointing initially in the z-direction is perpendicular to 
the group velocity. Specifically, for an initial state of the 
form l|14|l moving along the y-direction (ko x = 0), one 
finds in the limit dko y 3> 1 



(*p\x H (tm = 



2h 



1 — cos 



2ak 3 



(43) 



Thus, the amplitude of the zitterhewegung is again pro- 
portional to the wave length \o y — 2ix /ko y of the particle 
motion perpendicular to it. The frequency of the zitter- 
hewegung is given by fl — 2ak,Q y /h. Winkler et al. |45j 
have studied both theoretically and experimentally the 
magnitude of the Rashba spin orbit coupling in GaAs- 
based quantum well samples with heavy-hole densities 
of a few 10 14 m -2 and have found typical values for the 
characteristic length scale ma/H 2 of a few nanometers. 
Assuming a value of 2nm this corresponds to a coupling 
parameter of a — 0.3eVnm 3 , where we have used the 
heavy- hole mass m w 0.5too for GaAs [3!|- For a typical 
wave vector with ko y w O.lnm -1 this leads to frequencies 
Cl of order 10 11 Hz, an estimate which is of a similar order 
of magnitude as for the case of an n-doped quantum well. 



IV. CONCLUSIONS 



using the notations p± = p x ±ip y , <j± = <r x ±ia y , where p, 
a denote the hole momentum operator and Pauli matri- 
ces, respectively. These Pauli matrices operate on the to- 
tal angular momentum states with spin projection ±3/2 
along the growth direction; in this sense they represent a 
pseudospin degree of freedom rather than a genuine spin 
1/2. In the above equation, m is the heavy-hole mass, 
and a is Rashba spin-orbit coupling coefficient due to 
structure inversion asymmetry. The components of the 



We have investigated the notion of zitterhewegung in 
both n- and p-doped III-V zinc-blende semiconductor 
quantum wells, extending previous work on the two- 
dimensional electron gas |2|. zitterhewegung in the two- 
dimensional electron gas has been studied for the case 
of spin-orbit coupling of both the Rashba and Dressel- 
haus type. In the context of these investigations we have 
also discussed recent work by Zawadzki [4| . The crucial 
difference between zitterhewegung of free electrons and 
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electrons bound in the above semiconductor nanostruc- 
tures is the fact that in the latter systems the frequency 
of the oscillations is by orders of magnitude smaller while 
the amplitude is grossly larger. This circumstances make 
such systems favorable candidates for the experimental 
detection of zitterbewegung. 

The case of electron dynamics in a quantum wire is 
studied in great detail for various initial conditions. For 
this type of system we propose possible experiments for 
detecting the zitterbewegung of electronic wave packets. 
For an harmonic quantum well in the presence of Rashba 
spin-orbit coupling, the dynamical parameters can be 
tuned to a resonance condition where the amplitude of 
the zitterbewegung becomes maximal Q. This property 
should facilitate the experimental observation of this ef- 
fect. In addition to the orbital dynamics, we have also 
analyzed in detail the electron spin dynamics, which also 
show peculiarities at the resonance. 



Finally we have also discussed in detail the zitter- 
bewegung in the two-dimensional hole gas. Here we 
have considered spin-orbit coupling of heavy holes due 
to structure-inversion asymmetry. Concerning the fre- 
quency and amplitude of the zitterbewegung, similar re- 
sults are obtained as for the two-dimensional electron 
gas. 
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FIG. 1: Zitterbewegung of an electron in a harmonic quantum 
wire perpendicular to the wire direction. The electron is ini- 
tially injected into the lowest subband with its spin pointing 
along the 2-direction. The wave number ko y for the motion 
along the wire is ko y X = 5. The amplitude of the zitterbewe- 
gung is maximal at the resonance 2ako y = Hui (middle panel) . 
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FIG. 2: Amplitude of the zitterbewegung perpendicular to the 
wire direction as a function of fi/w = 2akoy/fkd for different 
values of the wave number ko y along the wire. Again the 
electron spin points initially in the ^-direction. 
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FIG. 3: The spin expectation values {tl> z \Bii(t)\ipz) as a func- 
tion of time at the resonance condition 2ako y = huj with 
fcoyA = 5 (cf. mid panel of Fig. 0. The time evolution of 
(V'zl'^f/WIV'z) (solid line) is governed by a single frequency as 
shown in Eq. 138H . while the other components (^z\o v H {t)\%l) z ) 
and {V , z|f|r(*)|V , z) (dashed lines) show beatings between two 
frequencies given by fi± . 



